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Abstract 

Recent CDF measurements of the inclusive cross section for a double parton scattering 
attach a great importance to any theoretical calculations of two-particle distribution func- 
tions. Using a parton interpretation of the leading logarithm diagrams of perturbative 
QCD theory, generalized Lipatov-Altarelli-Parisi-Dokshitzer equations for the two-parton 
distributions are re-obtained. The solutions of these equations are not at all the product 
of two single-parton distributions what is usually applied to the current analysis as ansatz. 
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The CDF Collaboration has recently measured a large number of double parton scatterings 
[1]. Thus new and complementary information on the structure of the proton can be obtained 
by identifying and analyzing events in which two parton-parton hard scatterings take place 
within one pp collision. This process, double parton scattering, provides information on both 
the spatial distribution of partons within the proton, and possible parton-parton correlations. 
Both the absolute rate for the double parton process and any dynamics that correlations may 
introduce are therefore of interest. The theoretical estimations of the effect under consideration 
have been done in a number of works [2, 3, 4, 5, 6, 7, 8]. 

For instance, the differential cross section for the four-jet process (due to the simultaneous 
interaction of two parton pairs) is given by [4, 5] 

da^y^ ^^^^ Dp{xi, X3) Dp{x2, X4), (1) 

daij stands for the two-jet cross section. The factor cxe// in the denominator represents the 
total inelastic cross section which is an estimate for the size of hadron, a^/f — TrVp. Based on 
the simple "hard sphere" model of proton structure, the expected value for (Tg// is 11 mb and is 
consistent with CDF measured value of (14.5±1.7l2!3) mb [1]. The two-parton distributions are 
supposed to be the product of two single-parton distributions times a momentum conserving 
phase space factor 

Dp{xi, Xj) = Dp{xi, Q'^) Dp{xj, Q"^) (1 - - Xj), (2) 

where Dp{xi,Q'^) are the single quark/gluon momentum distributions at the scale (deter- 
mined by a hard process) . 

The main purpose of this Letter is to analyze the status of the factorization ansatz (2) for 
many parton distributions in the perturbative QCD theory. 

Here one should note that the generalized Lipatov-Altarelli-Parisi-Dokshitzer equations for 
many parton distribution functions have been derived for the first time in refs. [9, 10] within the 
leading logarithm approximation of QCD using a method by Lipatov [11]. Under certain initial 
conditions these equations lead to the solutions, which are identical with the jet calculus rules 
proposed for multiparton fragmentation functions by Konishi-Ukawa-Veneziano [12]. Because 
of a very old affair it is necessary to recall some features of that investigation to be clear. 

In ref. [13] the structure functions of ep scattering and e'^e" annihilation were calculated in 
the leading logarithmic approximation for vector and pseudoscalar theories. Similar calculations 
in QCD were made in [14]. Lipatov shown [11] that the results of these calculations admit a 
simple interpretation in the framework of the parton model with a variable cut-off parameter 
A ~ (5^ with respect to the transverse momenta, and derived an equation for the scaling 
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violation of the parton distribution Dj {x, A) inside a dressed quark or gluon, which is in fact 
equivalent to the one proposed by Altarelli and Parisi [15] , within the difference that it was not 
applied by Lipatov to QCD. 

Calculating ladder graph contributions to the structure functions it is convenient to use 
Sudakov's technique [16]: 

2 1 

k^aq' + Pp + k±, q' ^ q p, q'^ = 0, pq' = -s, (3) 

S Zi 

q is the momentum of the virtual photon, p is the momentum {j? — 0) of the considered quark 
or gluon. In deep inelastic processes this constituent is one of the partons found at the reference 
virtualness —q^ = fi'^ in the target proton. In e^e~ annihilation it is that final quark or gluon 
which converts into the detected hadron. q is a light-like vector combined from p and q. k± is 
the component of k orthogonal to q and p. In the next step one takes the residues with respect 
to a in the Feynman integrals. Now each line r carries only (3r and k±r- And so one may 
introduce the amplitudes ^"(/Sr, k±r) for finding in the dressed constituent i n bare partons of 
given types and parameters k±_ri r — 1, n. Analogously ^"(/3r-, k±_r) are the probability 
amplitudes for finding a bare parton of i-type in a state of n dressed partons with Sudakov 
parameters k_i_r in the case of e^e~ annihilation. These amplitudes can be calculated as 
tree diagrams in an old-fashioned perturbation theory and obey the normalization condition 

1 = + E / n '^e{^)cfk^rm?nY.k^r)5{Y.^ - 1). (4) 

n=2-' r=l Pr 

where Zi is the wave function renormalization constant, to be interpreted as the probability 
that the constituent i remains in the bare state. 

An important feature of the leading logarithm diagrams is the ordering of the lines. The 
absence of interference-type diagrams allows a classical probability interpretation. The integra- 
tion region in k±_ essential in the leading logarithm approximation goes up to A = |A;_l| < A. 
The same cut-off is used for the ultraviolet-divergent Zi. The unrenormahzed amplitudes , 

are connected with the renormalized (cut-off independent) ones (/?", (^f via 

3 

*r = (5) 

Uj being the number of partons j in the state n. 

In the normalization condition (4) the left-hand side does not depend on A, while on the 
right-hand side the functions and the limits of integration over k±_ depend on A. Differen- 
tiating (4) with respect to A one obtains (see refs. [11, 9, 10] for details) 

= A^^r^ + a;.(^i), (6) 
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where 



(7) 



d\n{A/\k^\) 

Hi is the one-loop self-energy part of the constituent i, g\ is the running couphng constant. 
In eq. (7) one differentiates only the upper limit of the loop integral. 
The parton distribution function can be defined in this framework as 

DI{X,A) = ZAAX - 1) + E E / n -^0{/3r)d'k^r ■ 

n=2r{j) r=l P'r 

•l*ri'<^'(EM<^(E/5^ - - (8) 

The second sum X] runs over all partons of type j in the state n. The result of differentiation 
reads 

dDl{x,A) gl ^ fdx'f, , fx' 



r X 



where 



X 



E ^3-*hr + E ^3 — E ^3-^3132^ (10) 

jl<j' j'<jl jl<j2 

ujj^j^j,^{x — ^ Xi)dxidA/ A, as usual, is the probability that a parton of type j with the fraction x 
of the longitudinal momentum decays into two partons of types ji and 72 1 one of which has the 
fractions Xi of the longitudinal momentum and the transverse momentum A. These probabilities 
are defined as in eq.(7) but with the types of partons and the longitudinal momentum fractions 
in the loop specified. Also, the eq.(6) was used to derive the evolution equation (9). 
By introducing the natural variable 



1 , 
— -In 
27r6 



33 -2n/ . 
h = m QCD, 



one obtains 

dDl{x,t) ^ f dx' , fx 



i:l^-^Di\.'.t)p,j^]. (11) 



dt 

3 X 

This is an equation in the form of Altarelli and Parisi [15]. It is interesting that expression (10) 
for the kernels P already includes a regularization at x ^ x', which was introduced in ref. [15] 
afterwards. 

Defining the two-parton distribution function 

Df^'(xi,a;2, A) = ZiSij^Sjj^j^Sixi - X2)S{xi - 1) - 5ju2^(^i - a;2)-Df (xi. A) + 

+ EE E / n#W)rf'^±.l*rr<^'(EM<^(i-E/3^)W(.i)-^i)Wo.)-^2), (12) 

"=2 r(ji) r(j2) r=l Pr 



one obtains by differentiating with respect to A 

h' XI ^ ^ 



J2' X2 

where the kernel 

9^ 1 fxi\ 

■^-Pj-^hh \ — ]^ ujj^nj, (x xi) + u;j^nj2 (x xi) (14) 
oTT X \ / Ji^J2 n>n 

is defined without 5-function regularization. This is the generahzed Lipatov-AhareUi-Parisi- 
Dokshitzer equation for two parton distributions Dl'^-''^{xi,X2,t), representing the probabihty 
that in a dressed constituent i one finds two bare partons of types ji and j2 with the given 
momentum fractions xi and X2 ■ The result for the m-parton functions can be found in ref. [9]. 

It is readily verified by direct substitution that the solution of eq. (13) can be written via 
the convolution of single distributions [9, 10] 

Dr(xi,X2,i)= (15) 

This coincides with the jet calculus rules [12] proposed originally for the fragmentation func- 
tions and is the generalization of well-known Gribov-Lipatov relation installed for single func- 
tions [13, 14] (the distribution of bare partons inside a dressed constituent is identical to the 
distribution of dressed constituents in the fragmentation of a bare parton in the leading loga- 
rithm approximation) . The equations for the multiparton fragmentation functions are obtained 
by Lipatov's method in a similar way [10] and beyond the given investigation. 

The solution (15) shows that the distribution of partons is correlated in the leading logarithm 
approximation : 

D^^-{x^,X2.t) ^ D^{x,,t)Dr{x2.1^. (16) 

Of course, it is interesting to find out the phenomenological issue of the equations under con- 
sideration. It can be done within the well-known factorization of soft and hard stages (physics 
of short and long distances) [17]. As result the equations (11) and (13) describe the evolution 
of parton distributions in a hadron with t (Q^), if one replaces the index i by index h only. 

However, the initial conditions for new equations at t = (Q^ = /x^) are unknown a priori 
and must be introduced phcnomcnologically or must be extracted from experiments or some 
models dealing with physics of long distances (at the parton level: Dl{x, t = 0) = dijd{x — 1); 
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Dj^-'^{xi,X2,t = 0) = 0). Nevertheless the solution of eq. (13) with the given initial condition 
may be written as before via the convolution of single distributions [10] 



Di^^^ix„x2,t)= E /^T^i^P'(;^i,;^2,0)Dj^,(f^,t)z^jMt'^) + (17) 

J'31J2' a XI X2 \ J 

The reckoning for unsolved confinement problem (physics of long distances) is unknown 
two-parton correlation function D]^ ^'^ (2:1, 2^2, 0) at some scale //^. One can suppose that this 
function is the product of two single-parton distributions times a momentum conserving factor 
at this scale //^: 

^^(^1,^2,0) = D^^(^i,0)Di^(^2,0)^(l-^i-^2), (18) 

then 

L>r(xi,X2,t) = {Di^{x^,t)Di-{x2,t) + (19) 
E I^l'it D^^'(;^i,0)Di^'(;^2,0)DjJ,(f^,t)I^j,Mt,t)[^(l - ^1 - z,) - 1])9{1 - x^-x,) + 

t 1 l-xi ., / \ • 

. E. / rft' / ^ / ^ Di {z, + .2, t')j^^Pr^,,,, ^ Dj^,i^, t - t')D^l,{^, t-t% 

]'3i 02' Q XI X2 \ J 

where 

dz , . ,X 



Di{x,t)=j:f^D{{z,0)DjAt) (20) 



3' X 



is the solution of eq. (11) with the given initial condition Df^{x, 0) for parton distributions inside 
a hadron expressed via distributions at the parton level. 

This result (19) is, as matter of fact, the answer to the question set. If the two-parton distri- 
butions arc factorizcd at some scale fi^ then the evolution violates this factorization inevitably 
at any different scale (Q^ 7^ /i^), apart from the violation due to the kinematic correlations 
induced by the momentum conservation (given by ^^-functions)^. 

For a practical employment it is interesting to know the degree of this violation. It can be 
done numerically using, for instance, the CTEQ-parametrization [18] for single distributions 
as an input in eq. (19) and considering the kinematics of some specific process. Partly this 
problem was investigated theoretically in refs. [10, 19] and for the two-particle correlations of 
fragmentation functions in ref. [20]. That technique is based on the Mellin transformation of 
distribution functions like 

1 

Mi{n,t) = I dxx"Di{x,t). (21) 



^Here this is the analog of a momentum conserving phase space factor in eq. (2) 
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After that the integro-differential equations (11) and (13) become systems of ordinary hnear- 
differential equations of first order with constant coefficients and can be solved exphcitly [10, 19]. 
In order to obtain the distributions in ^-representation an inverse MeUin transformation must 
be performed 

Di{x,t) = J^.x--Miin,t), (22) 

where the integration runs along the imaginary axis to the right from all n-singularities. This 
can be done numerically again. However the asymptotic behaviour can be estimated. Namely, 
with the growth of t [Q"^) the second term in eq. (17) becomes dominant for the finite xi and 
X2 [19]. Thus the two-parton distribution functions "forget" the initial conditions unknown a 
priori and the correlations perturbatively calculated appear. 

The CDF Collaboration found no evidence for the kinematic correlation between the two 
scatterings in double parton events. This can mean only that the factorization ansatz (2) is 
the acceptable approximation at the scale Qcdf ~ 5 GeV accessible to CDF measurements 
{E?^^{min) ~ 5 GcV [1]). There are no arguments to assert that this ansatz (2) is acceptable 
at the larger scales of hard processes accessible to LHC measurements. The second term in 
eqs. (17) and (19) can give the meaning contribution to the cross section of many jet production 
at LHC energy and influence on the background estimations to the Higgs production. 

To summarize, the analysis shows that within the leading logarithm approximation of the 
perturbative QCD theory and the factorization of physics of short and long distances, the two- 
parton distribution functions being the product of two single distributions at some reference 
scale become to be dynamically correlated at any different scale of a hard process. These 
correlations are perturbatively calculable (19). Of course, in order to be more conclusive one 
needs to do numerical estimations of this effect that is planning to perform in the future. 
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